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We extend results developed by Chandler [J. Chem. Phys. 65, 2925 (1976)] for the dielectric 
constant of neutral site-site molecular models to mixtures of both charged and uncharged molecules. 
This provides a unified derivation connecting the Stillinger-Lovett moment conditions for ions to 
standard results for the dielectric constant for polar species and yields exact expressions for the 
small-fc expansion of the two-point intermolecular charge-density function used to determine the 
total Coulomb energy. The latter is useful in determining corrections to the thermodynamics of 
uniform site-site molecular models simulated with spherically truncated Coulomb interactions. 

In this note we extend results for the small wavevector expansion small-fc relationships between the dielectric 
constant and moments of the intermolecular and intramolecular correlations functions, originally developed by Chan- 
dler [l[ • This extension generalizes these relations to mixtures of neutral and charged molecular species and connects 
them to the small wavevector expansion of the two-point intermolecular charge-density function used to determine 
the total Coulomb energy. These results are used elsewhere to derive corrections to the thermodynamics of uniform 
site-site molecular models simulated with spherically truncated Coulomb interactions Q. 

I. DERIVATION OF ZEROTH AND SECOND MOMENT CONDITIONS FOR A MIXTURE OF 
NEUTRAL AND CHARGED SITE-SITE MOLECULES 

The total Coulomb energy obtained during simulation of a mixture of small site-site molecular species without any 
intramolecular charge-charge interactions is 

/ -, Nm n m> n M n M> \ 
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In this notation, the angular brackets indicate a normalized ensemble average, M and M' indicate a given molecular 
species, i and j indicate a given molecule of a given species, and a and 7 represent the intramolecular sites [l|, |3|. The 
Kronecker deltas are necessary to exclude any charge-charge interactions between intramolecular sites within a given 
molecule. This energy U q can be more compactly represented as 

where p qq is a two-point intermolecular charge-density function that explicity excludes any purely intramolecular 
charge correlations, as implied by Eq. (fTJ) and detailed below. 

The composite function p qq (r,r') is a charge-weighted linear combination of all intermolecular, two-point site-site 
distribution functions, 

p"(r,r') = 2J 51 q a Mq~tM> PaM~fM> (3) 
aM jM' 

We now relate this function to the basic charge-charge linear response function used in the theory of the dielectric 
constant. 
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For solutions of primitive model ions, charge neutrality and screening place specific requirements on the behavior of 
p qq in fc-space at small k values More generally, for a fluid composed of charged and polar molecules, the dielectric 
screening behavior of the molecules places restrictions on the decay of this two-point charge density. Based on this 
observation, we are able to harness a theoretical development of Chandler [lj that expresses the dielectric constant in 
terms of an exact sum of charge-density- weighted pair correlation functions. We generalize the derivation to include 
both charged and neutral site-site molecules and we take the dielectric constant as a given. From this vantage point, 
we may instead use these relations to place requirements on the decay of the two-point charge density p qq . 

We first define the instant ant eous single-point charge-density p q (r, R), a function of both a given external position 

r and the set of positions of all mobile particles R = {Rim - } = | r a/|? as 

N M ri M 
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With such a definition, the ensemble-averaged charge-density profile p q (r) is 

p q (r) = (p q (v,n)). (5) 

In the case of a uniform system, p q (r) — 0. Comparing Eqs. |T]) and |2]) and using Eq. ((4]), we may also express 
p qq (\r — r'|) for a uniform system as 

/ rtM n.M \ 
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We have used the equivalence of all molecules of type M in the last term. This term removes purely intramolec- 
ular charge-density correlations; we shall determine the small-A: contributions from that term based on well-known 
molecular properties using the approach of Chandler [l[ later in this note. The first term, in contrast, is exactly the 
charge-charge linear response function for a uniform neutral system: 

(p*(r,R)pV,H)) = (Sp q (r,n)Sp q (r',K)) = X qq (\r - r'|). (7) 

Here 8p q {r 1 R) = p q (r, R) — (p 9 (r, R)). Physically x qq describes the coupling between charge-density fluctuations at 
positions r and r'. As is well established |5|, such a function is intimately related to the dielectric behavior of the 
fluid at long distances, and furthermore, may be easily analyzed using basic electrostatics and standard definitions of 
the functional derivative. 

The electrostatic potential at r induced by a fixed external charge distribution Poxt( r ') ( e -ff-i a test charge Q placed 
at the origin, as considered by Chandler) is given by 

Vext(r) = j f^dr\ (8) 

and the associated electrostatic energy for a particular microscopic configuration characterized by the set of molecular 
positions R is then 



' p 9 (r,R)V oxt (r)dr. (9) 



This energy contribution will appear in the nonuniform system's Hamiltonian when V ox t(r) is nonzero. As such, we 
know from standard definitions of functional differentiation of free energies [5|, [6( that 

-<✓('.=)>»-/•«'>. a»> 

where j3 = (keT) -1 and the subscript V indicates that the ensemble average is taken in the presence of an external 
potential. Similarly we have 
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Xl q (r,r'). (11) 
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The total electrostatic potential at position r in the nonuniform fluid is then given by the sum of the external 
potential and the induced polarization potential: 

V to t(r) = Vext(r) + V P oi(r) = V cxt (r) + J dr' (12) 

To get a formula for the dielectric constant we expand about the uniform neutral system and evaluate V po i to linear 
order in V ox t using Eq. (jlip and find 

V tot (r) « V ext (r) + / ^ / dr" h^r")] , (13) 

= V cxt (r) - J dr' — J dr" f3 X qq (|r' - r"|) V cxt (r"). (14) 

Here x qq is the linear response function in the uniform fluid as in Eq. |7|. Taking the Fourier transform of the final 
equation, we find 

Vtot(k) = V cxt (k) - ^(3 x qq (k)V cxt (k). (15) 



Thus to linear order we have 

-X qq (k). (16) 



V to t(k) _ 1 47T/3 
Vext(k) k 2 

Phenomenologically, we know that in the limit of k — > 0, this ratio of the total electrostatic potential to the 
externally-imposed potential is exactly 1/e. Therefore, we find for our molecular mixture the general result 

te^-^wH- (17) 

Based on the limit in Eq. (|17[) . and expanding \ qq for small k as x (0)w + X (2)qq 

k , we have 

X {0)qq = 

4TrPx {2)qq = 1 - -• (18) 
e 

Any mixture with mobile ions acts as a conductor with e = oo in Eq. (|18p , independent of the nature of the neutral 
components. 

Our goal is to write a small-A: expansion of the two-point intermolecular charge density, 

p qq {k) « p ( ^ qq + k 2 p^ qq + 0(k 4 ). (19) 
As stated at the beginning of this derivation, 

p qq (\r - r'|) = X qq (\ r -r'\)-(j2N M J2J2 ^Ml^i* ^)5[v' - r« ) ) . 

\ M a=l7=l / 

We already have results for the small-A: expansion of the charge-charge linear response function \ qq ■ Now we must re- 
move the intramolecular contributions as described by the second term in the equation above. Defining the conditional 
singlet intramolecular site density functions £> a | 7 M( r l r ') for a / 7 as 

P 7 M(r'K| 7 A/(r|r') = (N M 5(r - rfi>)5(r> - r& )) , (20) 
and applying consequences of uniformity, we find 

p qq (\r - r'|) = X qq (\r - r'|) - J>m ]T QcMQ-yM [S ay S(r - r') + g ahM (\r - r'|)] . (21) 

M a, 7 
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Using the notation of Chandler for the term in square brackets [1| , we may write this expression more compactly as 

p«(|r - r'|) = x" (|r - r'|) - E feM<Z 7 M^ 7 M(|r - r'|). (22) 

M ct,7 

For neutral molecules, Chandler demonstrated that the small-A: components of Co ai M{k) are related to simple properties 
of the molecule. For both charged and uncharged molecules, the zeroth moment of Co ai u is simply 

^ajM = S <*y + J dr £ a | 7M (r) = 5 ai + (1 - 6 Q7 ) = 1. (23) 
Using this exact expression in Eq. ([22]) yields 

p i0)Qq = X (0)qq - E ^ E 9«M9 7 M = - E PM«M > (24) 
M Q,7 M 

an expression encompassing the standard zeroth moment condition for ions [4j and the zeroth moment for neutral 
molecular species [l|. 

(2) I — h 

The expression for 6j m determined by Chandler [1[ may be written most generally as 



to 



(2) 
M 



M 

1 f 1 

E QaMg-yMU^M = ~ g / dr E laMQ-yM Qa^M ( r ) r2 = _ g E ^aMQyM (^ 7 m) > ( 25 ) 

where Zq, 7 m is the bondlength between sites a and 7 for a molecule of species M. As shown in Ref. U for a neutral 
molecule indicated by iV below, the final summation in the above equation is simply related to the molecular dipole 
moment pn and the molecular polarizability as 

= E laNQiN^N = gMjV + k B Ta JV- (26) 

This relationship does not hold for a charged molecule since the dipole moment then depends on the choice of 
coordinate system. 

Distinguishing charged species (C) and neutral species (N) where {M} — {N}U {C}, and without substituting for 
toff and u>c , we find 

~(2) qq k B Te-l V"^ -(2) k B T e - 1 ^ .(2) v-^ -( 2 ) fOT\ 

P = ^ l^PM 2^ <laMq-yMUJ^ M = ^ P NUJ N ~ 2^ P cUJ C ■ ( 27 ) 

M a^j " N C 

~ (2) ~ (2) 

Thus, we may write a general expression for p qq in fc-space. Utilizing the expressions for ui N and lo c , we have 

P qq (k) = E PC& + k 2 ^—-k 2 Y.PN\\p 2 N + k B Ta w } +fc 2 ^E^E ^cq-yC <£ 7 c) + O (fc 4 ) . (28) 

Unlike x" m Eq. (JTHJ) , we see that the small-fc behavior of this two-point intermolecular charge-density function 
depends on several simple properties of the solution as a whole, like the dipole moment and polarizability of individual 
neutral molecules, and the net molecular charge and the average square bond lengths of charged molecules, as well 
as the dielectric constant. Thus by knowing simple single molecule properties and the long wavelength dielectric 
constant, we know how intermolecular charge-charge correlations decay in solution. This is the essential equation 
used to develop energy and pressure corrections for simulations of bulk liquids using molecular models with truncated 
Coulomb interactions 0. A related expression may be developed for larger molecular species with intramolecular 
charge-charge interactions. 
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